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We consider a ternary mixture of hard colloidal spheres, ideal polymer spheres, and rigid vanish- 
ingly thin needles, which model stretched polymers or colloidal rods. For this model we develop a 
geometry-based density functional theory, apply it to bulk fluid phases, and predict demixing phase 
behavior. In the case of no polymer-needle interactions, two-phase coexistence between colloid-rich 
and -poor phases is found. For hard needle-polymer interactions we predict rich phase diagrams, 
exhibiting three-phase coexistence, and reentrant demixing behavior. 



I. INTRODUCTION 

The richness of phase behavior of systems with purely 
repulsive interactions depends crucially on the number 
of components. For a one-component system like col- 
loidal hard spheres, there occurs a freezing transition 
from a single fluid phase to a dense crystal. Adding a sec- 
ond component, such as non- adsorbing globular polymer 
coils [Q, or rod- like particles ||, [sj generates an effective 
depletion-induced attraction between colloidal spheres, 
leading to the possibility of demixing. This transition 
is an analog of the vapor-liquid transition in simple flu- 
ids: The phase that is concentrated in one of the com- 
ponents corresponds to a liquid, while the dilute phase 
corresponds to a vapor, and one frequently refers to such 
phases as colloidal liquid and colloidal vapor, although 
the "vapor" is concentrated in the added component. 

Generic theoretical models for such systems are those 
introduced by Asakura and Oosawa (AO) and indepen- 
dently by VrijH, g], Bolhuis and Frenkel (BF)|], and 
Widom and Rowlinson (WR)[^. The AO model com- 
prises hard colloidal spheres mixed with polymer spheres 
that are ideal amongst themselves but cannot penetrate 
the colloids. The BF model adds stiff vanishingly thin 
needles to a hard sphere system. Because of their van- 
ishing thickness, the needles do not interact with one 
another. Clearly, both models are similar in spirit, as a 
non-interacting component is added to hard spheres. In 
the WR model this is different; two species of spheres 
interact symmetrically, such that hard core repulsion oc- 
curs only between particles of unlike species. Hence a 
pure system of either component is an ideal gas. All of 
these model binary mixtures exhibit liquid-vapor phase 
separation, well-established by computer simulations and 
theories @, |, |, 0, [l^, |l| . The WR modelg |l| 
has been studied with a range of approaches, including 
mean- field theory (MFT)l[l5), Percus-Yevick (PY) inte- 



gral equation theory p^, |16|, pTf, scaled-particle th eory 
(SPT) as well as computer simulations |o|. 
The precise location of the liquid- vapor critical point was 



located by simulations about 50 percent higher than pre- 
viously thought |l^ , still a challenge for theories (for 
a recent integral-equation closure, see Ref. [0). 

In the AO and BF cases a reservoir description has 
proven to be useful. The reservoir density of either poly- 
mers or needles rules the strength of effective attraction 
and hence plays a role similar to (inverse) temperature in 
simple substances. Although the WR model features an 
intrinsic symmetry which seems to preclude such a de- 
scription, an effective model can also be formulated []l5|. 
In the present work we consider the phase behavior of 
a mixture of spheres, polymers and needles, a natural 
combination of the above binary cases. We note that 
our ternary model may provide insight into certain real 
systems, such as paints, which contain colloidal latex and 
pigment particles, polymer thickeners and dispersants, as 
well as many other components pT[ . 

Density functional theory (DFT) ||2^ is a powerful ap- 
proach to equilibrium statistical systems, possibly under 
influence of an external potential. Building on Rosen- 
feld's workp3[, a geometry-based approach was recently 
proposed that also predicts bulk properties, without the 
need of any input, allowing the AO|Q, BF[^, and 
WRp^ models to be treated. Here we combine these 
tools to derive a DFT for ternary systems. 

In Sec. II we define the model ternary mixtures of 
spheres, polymers, and needles. In Sec. HI the DFT is 



developed. Application to bulk phases in Sec. IV yields 
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the phase behavior. We finish with concluding remarks 
in Sec. 



II. THE MODEL 

We consider a mixture of colloidal hard spheres 
(species C) of radius Rc, globular polymers (species P) 
of radius Rp, and vanishingly thin needles (species N) of 
length L, with respective number densities pc(r), pp(r) 
and pAr(r, r2), where r is the spatial coordinate and f2 
is a unit vector pointing along the needle axis (see Fig. 
|l|). The pair interaction between colloids is Vcc = oo 
if the separation r between sphere centers is less than 
2i?c and zero otherwise. The pair interactions between 
like particles of both other components vanish for all dis- 
tances: Vpp = Vnn = 0. For polymers this is an as- 
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sumption strictly valid only at the theta point; for nee- 
dles it becomes exact in the present limit of large aspect 
ratio, where overlapping needles contribute a negligible 
fraction of configurations. The colloidal spheres inter- 
act with both other components via excluded volume: 
The pair interaction between colloids and polymers is 
Vcp = oo if r < Rc + Rp, and zero otherwise; the in- 
teraction between colloids and needles is Vcn = cxd, if 
both overlap, and zero otherwise. What remains to be 
prescribed is the interaction between needles and poly- 
mers. We consider two cases: i) ideal interactions such 
that VpN — for all distances, and ii) excluded volume 
interactions such that VpN = oo if needle and polymer 
overlap, and zero otherwise. We denote the sphere di- 
ameters by ac ~ '2Rc, crp — 2Rp, the sphere packing 
fractions by rjc — AnRfjpc/S, t]p = AnRppp/S, and use 
a dimensionless needle density p*pf = pnL^ ■ 



thermore, for needles, we follow]^ to obtain 

<(r,f2) = - Al5{v + ill), (4) 

4 J-L/2 

w^{y,Q.) = i [5(r + flL/2) + (5(r - , (5) 

and r is the needle center of mass. The function Wi 
describes the linear extent of a needle, whereas Wq is 
characteristic of its endpoints. For vanishingly thin nee- 
dles, both surface and volume vanish, and so do the 
corresponding weights, W2 = = 0. Technically, 
the Mayer bond is generated through ~fgN{r,ft) = 
wf (r) * Wq (r, rj) -I- 102^ {r, n) * {r, il), where r is the 
difference vector between sphere and needle position. 



B. Weighted densities 



III. DENSITY FUNCTIONAL THEORY 
A. Weight functions 

We start with a geometrical representation of the 
particles in terms of weight functions u;*^, where p = 
3,2,1,0 corresponds to the particles' volume, sur- 
face, integral mean curvature and Euler characteristic, 
respectively [p7|, and i = C,P,N labels the species. 
We will use 5 as a unifying symbol for the spherical 
species C and P, and denote the radius as i?, where 
R = Rc,Rp for S = C,P, respectively. The weight 
functions are determined to give the hard core Mayer 
bonds fij — exp{Vij ) — 1 by a linear combination of terms 
wl^{r) * W3_^(r), where the star denotes the convolution, 
f;(r) * h{r) — J d^x g{x)h{r — x). 

For spheres, the usual weight functions p^. p8[ are 



'(r)= e{R-r), wlir) =d{R-r), 



(r) r/r. 



(r) 



'(r)[rr/r^ 



(1) 

i/3}2) 



The weight functions are used to smooth the possibly 
highly inhomogeneous density profiles by convolutions. 





= pc{r) * w^{r), 


(6) 




= pp(r) * (r). 


(7) 




= pcir)*w^''{r,n), 


(8) 




= pp(r)*wr^(r,f2). 


(9) 


n^(r,0) 


= /9Ar(r, r2) * (r, 


(10) 



where v = 0, 1, 2, 3, wl, v2, 77i2, and r — 0,1; pc{i^), 
Pp{y) and Pn{y,^) are the one-body density distribu- 
tions of spheres, polymers and needles, respectively. Note 
that n^,n^,n-^ are "pure" weighted densities, involving 
only variables of either species |2S 
and 



,PN 



In contrast, rig'^ 
are a convolution of the sphere densities with 
orientation-dependent weight function, combining char- 
acteristics of both species p5|. 



C. Free energy density 



where r — |r|, S{r) is the Dirac distribution, 
6{r) is the step function, and 1 is the identity 
matrix. Further linearly dependent weights are 
(r) = u;|(r)/(47ri?),w^i(r) = w^2(r)/(47ri?), w^(r) = 
wf (r)/i?. Note that these weights have different tenso- 



rial rank: 



^1 J 



^2 : 



are scalars; w^;^, are 



vectors; is a (traceless) matrix. These functions 

give the Mayer bond between pairs of spheres [Esl through 



-/ss/2 = * -h * wf 



w,^9 * w,^i . However, 



''v2 "vl- 

they are not sufficient to recover the sphere-needle Mayer 
bond[p7[. This is achieved through 



(r,0) = 2|w^2(r).f2 



(3) 



which contains information about both species: it is non- 
vanishing on the surface of a sphere with radius R, but 
this surface is "decorated" with an r2-dependence. Fur- 



The Helmholtz excess free energy is obtained by inte- 
grating over a free energy density, 

Fc.c[pc,Pp,Pn] - fceT J d^x J ^ $ ({n^}) , (11) 

where /cb is Boltzmann's constant, T is temperature, and 
the (local) reduced excess free energy density $ is a sim- 
ple function (not a functional) of the weighted densities 
This leads to a dependence of $ on orientation and 

, and n 



position. The variable x runs over space 1 23 
over the unit sphere |2^. 

The functional form of $ is obtained by consideration 
of the exact zero-dimensional excess free energy. We ob- 
tain 



$ = $ 



cc 



CP 



CN 



PN, 



(12) 
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where in the case of ideal polymer-needle interaction 
A = 0, and for hard polymer- needle interaction A = 1. 
In the following, the arguments of the weighted densi- 
ties are suppressed in the notation; sec Eqs. (|6|)-(|l0|) for 
the explicit dependence on r and fl. The hard sphere 
contribution, being equal to the pure HS case [Esl ESl, is 



73- 



-3detn^,2)/2 / [87r(l 



(13) 



The contribution due to interactions between colloids and 
polymers is the same as in the pure AO case|Q and is 
given by 



d<^cc p 



(14) 



The contribution due to interactions between colloids and 
needles is 



$CJV--n^ln(l-n^)-f 



-,N„CN 



1 



(15) 



Note that the simultaneous presence of <i>cp and ^cn 
in $ does not generate artificial interactions between P 
and N . For vanishing PN pair potential one can derive 
these terms from consideration of multi-cavity distribu- 
tions like in the binary CP (l|, ||] and CN cases H. In 
order to model the WR type interaction between poly- 
mers and needles in the presence of the colloidal spheres 
we use 



PN 



1-n? 



(16) 



This can be derived as follows. The starting point 
is a functional for binary hard spheres with added 
needles. Linearization in one of the sphere densities 
(which becomes the polymer species) is performed in the 
same way as linearization of binary hard spheres leads 
to the CP functional |]l^. In the absence of colloids, 
we obtain <i> = <i>cAr = 11^71^ + nfn2^. Then the 
density functional then can be rewritten as F^xc — 
-Jd^rJ d^r' J d'npp{r)fpNir; r', f2)pAr(r, J^)/(47r). 
This is precisely (a generalization to needles of) the 
mean-field DFT for the WR model Ahhough this 
does not feature the exact Od limit, as the geometry- 
based DFT^^ for WR spheres does, we expect 
differences to be small. 



IV. RESULTS 



A. Bulk fluid phases 



The hard sphere contribution is equal to the Percus- 
Yevick compressibility (and scaled-particle) result, which 



cc 



377c[3??c(2 - 7jc) - 2(1 - r^c? ln(l - yc)] 



(17) 

The colloid-polymer contribution is equal to that pre- 
dicted by free volume theory and rederived by 
DFTil as 



CP 



r/p/(8^i?3) 



{3q?7c [6(1 - vc)^ + 3q(2 -rjc- vl) 



+2q2(l + nc+ Til)] - 6(1 - T^c f ln(l - vc)) , (18) 

where q = ap/ac- The colloid-needle contribution 
equals the perturbative (around a pure hard sphere fluid) 
treatment of Ref. |^ , which can be shown to equal the re- 
sult from application of scaled-particle theory [^, and 
DFT pi, and is given by 



CN 



PN 



ln(l - j]c) + 



3L 



Vc 



ARc 1 - rjc 
The WR type polymer-needle contribution is 



PN 



1 



3£ \ PnVp 
ARp J 1 



VC 



(19) 



(20) 



For completeness, the ideal free energy contribution is 



$id = 



i=C,P.N 



p,[ln{p,Af) - 1], 



(21) 



where the A^ are (irrelevant) thermal wavelengths of 
species i. This puts us into a position to obtain the re- 
duced total free energy per volume <I>tot = 'I'id -I- $ of any 
given fluid state characterized by the bulk densities and 
relative sizes of the three components. 



B. Phase diagram 

The general conditions for phase coexistence are equal- 
ity of the total pressures ptot , and of the chemical poten- 
tials fj.i in the coexisting phases. Equality of temperature 
is trivial in hard-body systems. For phase equilibrium 
between phases I and II, 



Ptot 

pI 



-Ftot 



i = C,P, N. 



(22) 
(23) 



These are four equations for six unknowns (two state- 
points each characterized by three densities) . Hence two- 
phase coexistence regions depend parametrically on two 
free parameters. For three-phase equilibrium between 
phases I, II, and III 



For homogeneous density profiles, pi = const, the inte- 
grations in Eqs. (|6|)-(p^) can be carried out explicitly. 



Ptot 



-Ptot 



-Ftot 

,,in 



C^PN. 



(24) 
(25) 
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Eight equations for nine variables leave one free parame- 
ter. 

In our case ptot/ksT = -^tot + J2i=c,p,N Pid^tot/dpi, 
and Hi = kBTd^tot/dpi yield analytical expressions. We 
solve the resulting sets of equations numerically, which is 
st r aight forward . 



1. Ideal polymer-needle interaction 

Let us first explain our representation of the ternary 
phase diagrams. We take the system densities r]c,r]p, p*^ 
as basic variables. For given particle sizes, these span 
a three-dimensional (3d) phase space. Each point in 
this space corresponds to a possible bulk state, at some 
pressure ptot- Two-phase coexistence is indicated by a 
pair of points that are joined by a straight tie line. Ac- 
cordingly, three phase coexistence is a triplet of points, 
defining a triangle. In order to graphically represent the 
phase diagram, we show surfaces defined by one ther- 
modynamic parameter being constant. Such surfaces are 
conveniently taken such that coexistence lines (and trian- 
gles) lie completely within the surface. Clearly, this can 
be accomodated by imposing a constant value of ptot or 
any of fiCi P-Pi a-nd iin- Here we choose fip — const, and 
hence imagine controlling the system directly with rjc 
and p^, but indirectly via coupling to a polymer reservoir 
of packing fraction 77p = {An /3){Rp / Ap)^ exp{iip / ksT) . 
A constant-77p-surface is non-trivially embedded in the 
3d phase diagram. To depict it graphically, we show 
projections onto the three sides of the coordinate sys- 
tem, namely the rjc — rjc — rip, and rjp — p]^ planes, 
as well as a perspective 3d view. Furthermore, we in- 
dicate the accessible regions that are compatible with 
the constraint of fixed rfp. Their boundaries are implic- 
itly defined through rfp{ric — Q,rjp,p*j^) = const and 
VpivcVPT P*N = 0) = const. Note that tiehnes are al- 
lowed to cross inaccessible regions. 

For simplicity, and to establish a reference case, we 
initially ignore polymer-needle interactions and consider 
equal particle sizes, ac = crp = L. In the absence of 
polymer {rjp = 0), colloids and needles demix, as shown 
in Fig. |2^. Increasing the packing fraction of polymers 
in the reservoir causes the demixed region to grow and 
to shift to smaller rjc and (see Fig. ||b for rjp — 0.5). 
This behavior can be understood if addition of a second 
depleting species simply enhances the depletion-induced 
attraction between colloids. Increasing rjp further causes 
the critical point to hit the = axis. This is pre- 
cisely the demixing critical point of the binary CP (AO) 
model, which is located at rjp — 0.63831 (see Fig. ^). 
Computer simulations are currently being carried out to 
test the accuracy of this value For still larger rjp, 
the mixed states become disconnected, hence there is no 
path between colloid-rich and colloid-poor phases that 
does not pass through a first-order phase transition (see 
Fig. |d for 77^ = 0.8). 



2. Hard polymer-needle interaction 

Turning on the excluded volume interaction between 
polymers and needles allows the possibility of demixing 
between these components. In the absence of colloids, the 
PN mixture is of WR type: Interactions between parti- 
cles of like species vanish, while unlike particles interact 
with a hard core repulsion. Our case is a generalization 
to non-spherical particle shapes. In the mean-field treat- 
ment this does not affect the phase diagram, as only the 
net excluded volume enters into the theory. This robust- 
ness is also present in our approach. 

We first consider equal particle sizes, crc = up = L. 
It turns out that interesting behavior is observed only 
for small p^ < 0.2. The colloid-needle demixing curve 
lies well above this region, and is only weakly affected by 
rjp > 0. In the absence of needles (p^ — 0) and for large 
enough polymer density, colloids and polymers demix, 
indicated by a miscibility gap along the p^ = axis (see 
Fig. ||a for rjp = 0.8). Increasing needle density p^ > 
causes the gap to shrink and eventually to disappear in 
a critical point. Quite surprisingly, and in contrast to 
the former case of absent PA'' interactions, the addition 
of needles favors mixing. This behavior may reflect a 
competition between the depleting effects of interacting 
polymers and needles. By analogy with the CP subsys- 
tem it is clear that at sufficiently high polymer density, a 
PN miscibility gap will open for rjc = 0. However, this 
happens not by growing a small bump as in the CP case. 
Instead the CP demixing curve bends over to smaller rjc 
and touches (with its critical point) the rjc = axis (see 
Fig. ||b for rjp — 1.08731). For larger rjp, the critical 
point disappears (see Fig. ^ for rjp = 1.2). 

In order to bring CP and CN demixing closer together, 
we consider a reduced polymer size crp = ac/'^, generat- 
ing a weaker depletion attraction between colloids (at the 
same number density of polymers), and longer needles, 
L — 2(7 c generating stronger depletion between colloids, 
and hence lower ptf at the critical point in the binary 
CN case. Figure ^ shows the binodals in the (three) bi- 
nary subsystems. For the ternary mixture, we follow a 
path of increasing rjp, starting with rjp = 0, for which 
the phase diagram is displayed in Fig. There is no 
polymer present in the system, and phase separation into 
colloid-rich and needle-rich phases occurs at high enough 
densities of these components. Both jyp — p^ and rjc — rjp 
planes are inaccessible as rjp — 0. Increasing polymer 
density {rjp — 0.4 in Fig.||b) shifts the CN critical point 
to lower rjc, distorting the formerly rounded shape of 
the binodal. For rjc = 0, polymers and needles demix, as 
rjp is above the critical value for the Widom-Rowlinson 
type demixing of these species. The presence of colloids 
{rjc > 0) disturbs the PA-transition; the miscibility gap 
narrows, eventually disappearing in a critical point, with 
subsequent miscibility. At rj'p = 0.408107 (Fig. |c) the 
CN and PA critical points merge into a single one, and 
a needle-rich phase (A) becomes isolated. This coex- 
ists with a phase that consists (primarily) of colloids and 
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polymers at varying composition. For growing r]p, the 
"double" critical point broadens into a line and results 
in a thin neck joining both transitions. 

With increasing r]p the coexistence region broadens 
further (see Fig. |p for r]p — 0.5). Colloids and poly- 
mers also demix. For = 0, the system is above the 
critical point for the pure AO model, and hence coex- 
istence between colloid-rich and polymer-rich phases oc- 
curs. Again, the presence of the third component, in this 
case N, causes the density gap to shrink and eventually 
disappear with a critical point. As all binary subsystems 
are by now demixed, it is evident that the system will ul- 
timately display coexistence between three phases, each 
one enriched by one of the components, and represented 
by a triangle in system representation. Each corner of 
the triangle corresponds to one of the three coexisting 
phases. The Gibbs phase rule dictates that one degree of 
freedom remains, which is ptot or, equivalently, r]p (note 
that for hard-body systems, temperature is trivially re- 
lated to pressure). It is striking, however, how this tri- 
angle develops. One might expect this to occur by the 
joining of existing binary coexistence regions. This is not 
the case. The ternary region instead grows solely out of 
the A^-rich-poor coexistence, whereby CP-coexistence is 
only a spectator, separated by mixed states. The ini- 
tial three-phase triangle is extremely elongated (being a 
line as a boundary case). One corner corresponds to a 
needle-rich phase; both others differ only slightly in den- 
sities, one phase favoring colloids, the other polymers. 
Moving away from this CP-edge of the triangle (by re- 
ducing p^) leads to binary coexistence between C and 
P. This phase separation is reminiscent of the behavior 
of the pure AO model. However this reentrant coexis- 
tence is triggered by the presence of the needles, and it is 
separated (by mixed states) from the pure AO transition 
(and its region of stability in the presence of needles). In 
Fig. 1^ we show results for rjp — 0.52626, where the crit- 
ical points of both CP transitions have already merged, 
and again a neck is reminiscent of the formerly distinct 
transitions. For still larger r/p, the three-phase triangle 
grows further (see Fig. ^ for rjp = 0.54). Ultimately, at 
sufficient concentration the colloids must freeze, but we 
disregard the solid phase in the present work. We finally 



note that the whole scenario is covered over a relatively 
small density interval rjp = 0.4 — 0.54, and that the pack- 
ing fractions of colloids and polymers are only moderate. 
However, needle densities can be quite high. 



V. DISCUSSION 

In conclusion, we have considered a simple hard-body 
model for a mixture of spherical colloidal particles, glob- 
ular polymer coils and needle-shaped objects, which may 
represent either colloidal needles, stretched polymers or 
polyelectrolytes. We have extended a recent DFT ap- 
proach to this model and applied it to bulk fluid phases. 
The resulting phase behavior is very rich, ensuing from 
competition of demixing in the binary subsystems. 

The present work has interesting implications for the 
techniques of integrating out degrees of freedom (see e.g. 

^2[). Note that by integrating out, e.g., the needles, 
effective interactions between pairs of colloids, pairs of 
polymers, as well as colloids and polymers arise. Hence 
one arives at a binary mixture with (soft) depletion in- 
teractions. To what extent the ultimate mapping onto 
a one-component (colloid) system, by further integrating 
out the polymers, can be achieved is an interesting ques- 
tion. As a further outlook, the inclusion of freezing of 
colloids, disregarded in the present work, would further 
enrich the phase behavior. Computer simulations are de- 
sirable to test the theoretical phase diagrams. Further- 
more it is interesting to elucidate the structural correla- 
tions present in the various fluid phases. Inhomogeneous 
situations, such as induced by walls or present at inter- 
faces between demixed states, constitute further exciting 
directions of research. 
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FIG. 2: Demixing phase diagram of a ternary colloid- 
polymer-needle mixture with ideal polymer-needle interac- 
tions for ac = Op = L, and r/p =0 (a), 0.5 (b), 0.63831 
(c), 0.8 (d). Shown are binodals (lines), tielines between co- 
existing phases (thin lines), and critical points (dots). 




FIG. 3: Demixing phase diagram of a ternary colloid-polymer- 
needle mixture with hard polymer-needle interactions for 
ac = (TP = L, and =0.8 (a), 1.08731 (b), 1.2 (c). 



8 




9 




FIG. 5: Same as Fig. 0, but for oc ~ '2op — L/2, and rfp —Q 
(a), 0.4 (b), 0.408107 (c), 0.5 (d), 0.52626 (e), 0.54 (f). 



